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ABSTRACT  OF  THE  DISSERTATION 

Wave-Mean  Flow  Interaction  in  the  Storm-Time  Thermosphere 
Using  a  Two-Dimensional  Model 


by 


Douglas  George  Brinkman 
Doctor  of  Philosophy  in  Atmospheric  Sciences 
University  of  California,  Los  Angeles,  1990 
Professor  S.  V.  Venkateswaran,  Chair 


A  two-dimensional  pole-to-pole  numerical  model  with  background  sol¬ 
stitial  winds  has  been  used  to  study  the  global  dynamical  response  of 
the  thermosphere  to  high-latitude  energy  inputs  associated  with  a  model 
geomagnetic  storm.  This  model  storm,  consisting  of  a  series  of  sub¬ 
storms,  has  four  distinct  pulses  of  heat  input  over  a  12-hr  period. 

The  storm  generated  six  groups  of  wave  disturbances  in  the  form  of 
internal  gravity  waves.  Characteristics  of  these  disturbances  such  as 
their  vertical  structure,  phawe  propagation  speeds,  ard  influences  of 
dispersion  and  dissipation  are  illustrated.  The  disturbances  in  the 
velocity  field  and  in  the  temperature  field  are  shown,  along  with  their 
mutual  phase  relationship  with  regard  to  the  meridional  and  vertical 


transport  of  heat  and  momentum. 


In  addition  to  the  generation  of  disturbances,  the  thermospheric 
response  to  the  sustained  part  of  the  heat  input  is  \  to  consist 

of  the  establishment  of  a  global  meridional  circulation  which  is  ini¬ 
tiated  in  about  3  to  4  hours  after  storm  commencement  and  never  quite 
reaches  steady-state  in  the  simulation.  The  time  development  of  this 
circulation  and  its  associated  temperature  field  are  described.  - - *n 

The  main  purpose  of  this  study  is  to  investigate  the  interaction, 
between  the  disturbances  and  the  mean  flow  associated  with  the  storm. 

It  is  shown  that  this  interaction  can  be  represented  in  terms  of  an 
induced  circulation.  This  induced  circulation  is  forced  by  the  time- 
dependent  nature  of  the  eddy  flux  convergences  (divergences)  of  heat 
and  momentum.  The  influence  of  this  induced  circulation  is  detectable 
most  noticeably  in  the  heat  sources  and  sinks  set  up  by  the  Eulerian 
mean  circulation.  At  altitudes  above  150  km,  in  the  equatorial  region, 
the  equivalent  temperature  changes  due  to  such  sources  and  sinks  re¬ 
sulting  from  the  induced  circulation  are  1/3  to  1/4  of  the  changes  due 
to  the  Eulerian  mean  circulation. 

Spatial  and  temporal  variations  of  the  storm-time  winds  are  respon¬ 
sible  for  the  differences  between  Lagrangian  and  mean  Eulerian  trajec¬ 
tories.  Trajectory  calculations  also,  show  that  fluid  material  trans¬ 
ported  from  the  source  regions  does  not  reach  equatorial  latitudes. 
However,  strom-generated  waves,  reaching  the  equator  within  three  hours 
of  storm  onset,  initiate  fluid  motions  at  low-latitudes 


CHAPTER  I 


INTRODUCTION 


The  thermosphere  is  a  region  of  low  air  density,  end  high  values 
of  thermal  conductivity,  molecular  viscosity,  and  static  stability 
(Figures  1.1a  and  1.1b).  In  particular,  it  is  a  region  characterized 
by  its  quick  dynamic  response  to  energy  inputs.  A  major  source  of  en¬ 
ergy  for  the  thermosphere  is  solar  radiation  in  EUV  and  UV  wavelengths, 
which  is  variable  with  solar  cycle,  season,  and  time  of  day.  In  ad¬ 
dition,  there  is  an  important  source  of  energy  at  high  latitudes  re¬ 
sulting  from  the  dissipation  of  electrical  currents  in  the  auroral  re¬ 
gion  and,  to  a  lesser  extent,  from  the  precipitation  of  energetic  par¬ 
ticles  (Banks,  1977) .  This  high-latitude  energy  source  is  believed  to 
be  present  at  all  times,  but  its  strength  increases  considerably  dur¬ 
ing  geomagnetic  storms.  (Killeen,  1987;  Richmond,  1983;  Roble,1977.) 

The  global  dynamic  response  of  the  thermosphere  during  geomagnetic 
storms  has  been  the  subject  of  considerable  interest.  Earlier  studies 
have  indicated  that  this  response  consists  of  several  parts:  1)  A 
localized  zonal  jet  with  associated  vertical  circulations  in  the  region 
of  intensified  auroral  currents.  The  zonal  jet  arises  as  a  result 
of  collisions  of  electrically  neutral  air  molecules  with  ions  which 
are,  in  turn,  accelerated  by  polar  cap  electric  fields.  Such  electric 
fields  are  generated  by  magnetospheric  convection  and  mapped  onto  the 
ionosphere  along  highly  conducting  geomagnetic  field  lines.  Killeen 
et.  al.,  1984;  Killeen  et.  al.,  1986;  Larsen  and  Mikkelsen,  1983; 


1 


Mikkelsen  at.  al.,  1981;  Walterscheid  and  Boucher,  1984).  2)  Internal 

gravity  waves,  generated  especially  in  the  initial  phases  of  a  storm 
and  propagating  to  low  latitudes  (Francis,  1975;  Georges,  1968; 
Hernandez  and  Roble,  1978;  Hines,  1960,  1974;  Hunsucker,  1982; 

Richmond,  1978,  1979a;  Rees  et.  al.,  1984;  Roble  et.  al.,  1978;  Testud, 
1970)  .  3)  A  global  meridional  circulation  driven  by  the  sustained  heat¬ 
ing  maintained  at  high  latitudes  (Babcock  and  Evans,  1979;  Buonsanto 
et .  al.,  1989;  Burnside  et.  al.  1983;  Crowley  et.  al.,  1989b;  Hays  and 
Roble,  1971;  Hernandez  and  Roble,  1976,  1984;  Mazaudier  et.  al.,  1985). 

These  dynamical  responses  have  been  studied  using  both  two-  and 
three-dimensional  models.  Dickinson  et .  al.  (1975,1977),  and  Roble  et . 
al.  (1977),  used  a  time-dependent  two-dimensional  model  to  obtain  the 
steady-state  circulation  forced  by  solar  heating  and  high-latitude 
heating  for  various  levels  of  magnetic  activity.  They  also  investi¬ 
gated  the  effects  of  seasonal  and  solar  cycle  variations  in  maintaining 
the  steady-state  circulation. 

Richmond  and  Matsushita  (1975)  used  a  time-dependent  two-dimension¬ 
al  single-hemisphere  model  to  examine  the  response  to  an  isolated  sub¬ 
storm.  This  response  consisted  primarily  of  equatorward  propagating 
gravity  waves  along  with  strong  vertical  and  zonal  winds  at  high  lati¬ 
tudes.  In  later  papers,  Richmond  (1978,1979a)  used  this  model  to 
study  gravity  wave  generation,  propagation,  and  dissipation.  Further¬ 
more,  Richmond  (1979b)  made  a  quantitative  comparison  of  heating  due 
to  viscous  dissipation  of  gravity  waves  with  heating  due  to  the  meri¬ 
dional  circulation. 
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There  is  increasing  recognition  in  recent  years  that  the  full  com¬ 
plexity  of  the  storm-time  response  can  only  be  handled  by  three-dimen- 
ional  models.  Several  studies  using  models  such  as  the  NCAR-TIGCM  and 
UCL-TGCM  have  become  available.  (Crowley  et .  al.,  1989a;  Dickinson  et . 
al.,  1980;  Dickinson  et .  al.,1984;  Fuller-Rowell  and  Rees,  1981; 
Fuller-Rowell  and  Rees,  1984;  Roble  et.  al.,  1982;  Roble  et.  al.,  1983; 
and  Roble  et.  al.,  1987;)  However,  the  results  from  such  models  are 
complex  and  are  not  readily  interpreted.  Consequently,  two-dimensional 
models  may  still  play  a  complementary  role  in  the  diagnostics  of  specific 
results  from  three-dimensional  models.  Alternatively,  they  may  serve  as 
test-beds  for  preliminary  simulation  of  physical  processes  that  await 
inclusion  into  three-dimensional  models. 

The  present  paper  is  concerned  with  the  nature  and  importance  of 
wave-mean  flow  interaction  under  storm-time  conditions,  and  the  conse¬ 
quences  of  such  interaction  for  the  Lagrangian  and  Eulerian  descrip¬ 
tions  of  the  storm-time  circulation.  The  model  used  is  a  generaliza¬ 
tion  of  the  model  of  Richmond  and  Matsushita  (1975) ,  referred  to  ear¬ 
lier.  The  generalization  consists  of  the  extension  of  the  model  domain 
from  pole-to-pole  and  the  inclusion  of  seasonal  background  winds.  As 
a  consequence,  the  model  enables  us  to  study  solstitial  conditions 
when  the  winds  are  hemispherically  asymmetric.  The  influence  of  such 
winds  on  the  propagation  of  storm-induced  wave  disturbances  is  an  ad¬ 
ditional  interest  of  this  study. 

This  study  of  the  wave-mean  flow  interaction  for  the  thermosphere 
in  general  and  under  storm-time  conditions,  in  particular,  is  believed 
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to  be  new.  Ic  shows  that,  although  the  disturbances  generated  during 
the  storm  are  of  large  amplitude,  the  phase  speeds  of  the  disturbances 
(-500  m/s)  are  also  large  compared  to  the  background  wind  speeds  (50 
m/s)  so  that  critical  levels  (where  the  phase  speeds  match  the  back¬ 
ground  wind  speeds)  are  not  encountered  in  the  thermosphere.  This  is 
in  contrast  to  the  situation  at  lower  levels  in  the  atmosphere  where 
most  wave-mean  flow  interactions  occur  at  such  locations  (Booker  and 
Bretherton,  1967) .  However,  other  favorable  conditions  for  the  inter¬ 
actions,  like  wave  transience  and  dissipative  effects  (which  increase 
rapidly  with  height),  are  available  in  the  thermosphere.  The  strength 
of  the  interactions  is  ,  nevertheless,  moderated  by  the  large  vertical 
and  horizontal  scales  of  the  disturbances. 

The  layout  of  the  rest  of  the  dissertation  is  as  follows:  In  Chap¬ 
ter  2  we  give  model  details,  closely  following  the  work  of  Richmond 
and  Matsushita.  Chapter  3  presents  several  model  results  of  the  global 
circulation  and  of  the  waves  resulting  from  our  storm  model.  Chapter  4 
deals  with  wave-mean  flow  interactions.  Chapter  5  draws  a  distinction 
between  the  Lagrangian  and  Eulerian  mean  circulations  and  gives  exam¬ 
ples  of  Stokes  drift  calculations  for  several  locations.  Finally, 
Chapter  6  summeriz es  the  major  conclusions  of  this  study. 
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FIGURE  CAPTIONS  FOR  CHAPTER  I 


FIG.  1.1.  Typical  thermospheric  parameters:  a)  Temperature  (*K) 
{solid  line)  and  log  density  (KG/M^)  (dashed  line) .  b)  Coefficient  of 
molecular  viscosity  (X10~5  KG/MS)  (solid  line)  and  coefficient  of  mole¬ 
cular  thermal  conductivity  (X10“2  W/M’K)  (dashed  line) .  All  from 
Richmond  and  Matsushita  (1975),  Table  1. 
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CHAPTER  II 


TWO-DIMENSIONAL  THERMOSPHERIC  MODEL 

2.1  MODEL  EQUATIONS 

As  in  Richmond  and  Matsushita  (1975),  hereafter  abbreviated  as  RM 
our  model  equations  written  in  the  log-pressure  coordinate  system  are 
Zonal  momentum  equation: 

— 1 -  =  -  f  2ffcos9  +  -x-  cotG  Vo  +  X  (2.1) 

Dt  ^  a  J  ° 

Meridional  momentum  equation: 

DVn  1  J*  /  Va  n 

—r  m  "  A  I?  +  l  2Cfcos0  +  —  cote  J  Va  +  Y  (2.2) 

Dt  a  o0  ^  a  J  v 

Hydrostatic  equation: 


30 

3z 


=  RT 


(2.3) 


Continuity  equation: 


le  ( v0sin9 )  +  p  h  ( p” )  ■ 0 


Thermodynamic  energy  equation: 


D_ 

Dt 


(CpT)  +  RwT  »  Q  +  K 


(2.5) 
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and  the  Equation  of  State: 


pgH 


(2.6) 


The  meanings  of  symbols  in  the  model  equations  are: 
a  =  radius  of  earth, 

g  =  gravity  9.5  m/s^. 


P 

Ps 

t 

w 

2 

CP 

H 

R 

T 

<V*,V0) 

W 

z 

p 

«M) 

Q 


pressure . 

reference  pressure  at  80  km  (1.031Pa). 
time 

~  ,  a  measure  of  vertical  velocity,  positive  upward, 

height . 

specific  heat  at  constant  pressure, 
scale-height  (RT/g) . 
gas  constant. 


temperature . 

(zonal,  meridional)  velocity, 
vertical  velocity  (in  m/s) 
log-pressure  coordinate,  ln(ps/p). 
density. 

(longitude,  colatitude) . 


angular  velocity  of  the  earth. 


ao 

dZ 


RT. 


geopotential  (gz)  and 


D_ 

Dt 


+  + 
3t  a  30 


w 


3_ 

3z 


\r  addition,  X  and  Y  in  the  momentum  equations  are: 


1  3  (  \ 

"  P  3z  i  gT<*>  J  J<1 


and 


“  p  h  (  )  +  Je 


where 


(  ^ 

r  Mm  +  Ht  ^  3  f  mm  mm  } 

l  '  Te  J 

II 

l  H  3  92  V-  V4*  '  0  J 

and 

(  X(j)/  Iq  )  =  Ion  drag  in  the  (  4>,  0  )  directions* 

In  the  above,  and  (lt  are  coefficients  of  molecular  and  eddy  viscos¬ 
ity,  respectively.  A  typical  value  of  (1^  +  |lt  at  250  km  (middle  ther¬ 
mosphere)  is  4.0X10"*5  kg/ (m  s)  .  Values  at  other  heights  are  obtained 
from  Table  1  of  RM  (p.2841) . 

In  calculating  the  ion  drag,  the  model  assumes  that  the  ions  are 
not  influenced  by  electric  fields  of  magnetospheric  or  ionospheric 
origin.  They  do,  however,  collide  with  neutrals  and  provide  a  retard- 
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ing  force  on  air  motion. 


Similarly,  in  the  thermodynamic  energy  equation: 


and 


with 


Q 


K 


q 


where  Q  is  the  diabatic  heating  rate  (per  unit  mass)  split  into  three 
parts:  global  mean  (Qq)  due  to  solar  heating  only;  departure  from  the 
global  mean  solar  heating  plus  heating  from  high-latitude  auroral 
sources  (Qg ) ;  and  storm-time  auroral  heating  (QA)  .  K  is  the  heating 
rate  due  to  thermal  conductivity  consisting  of  both  molecular  (coeffi¬ 
cient,  Km  )  and  eddy  (coefficient,  Kt  )  conductivities.  A  typical  value 
for  molecular  conductivity  at  250  km  is  .09  W/ (m  *K)  .  Eddy  conductiv¬ 
ity  is  negligible  above  about  125  km.  Values  for  coefficients  of  mole¬ 
cular  and  eddy  thermal  conduction  are  shown  in  Table  1  of  RM  (p.2841) . 


2.2  MODEL  DESCRIPTION  AND  BOUNDARY  CONDITIONS 

The  model  is  two-dimensional  in  altitude  and  colatitude  and  extends 


from  80  km  (Z=0)  to  450  km  (Z=14.4)  in  the  vertical,  and  from  the  north 
(0=0*)  to  the  south  (0=180*)  poles  in  the  horizontal.  The  grid  spacing 


is  uniform  in  the  vertical  with  AZ=0.3. 


In  the  horizontal  the  coiati— 


tude  spacing  varies  from  1.2  degrees  near  the  poles  (for  a  more  detail- 


ed  representation  of  the  region  of  storm-time  heating)  to  2.4*  colati¬ 
tude  near  the  equator. 


The  boundary  conditions  for  the  model  are: 

V(J>  =  V0  -  0  at  0  =  0*  and  180*  . 

w  =  0  at  Z  =  14.4. 

z  -  80  km  at  Z  =  0 . 

T  =  180*  K  at  Z  =  0. 

The  only  difference  in  the  boundary  conditions  of  the  model  from  RM  is 
the  extension  of  the  model  domain  from  pole-to-pole . 


2.3  MODEL  INITIALIZATION 

The  model  thermosphere  is  initially  at  rest  and  has  the  global  mean 
temperature  field  shown  in  Figure  2.1.  As  in  RM,  this  temperature 
field  is  prescribed  according  to  Jacchia's  (1971)  model  atmosphere  with 
an  exospheric  temperature  of  1000*  K.  The  prescription  for  the  neutral 
atmospheric  composition  and  for  the  vertical  distribution  of  electron 
density  is  the  same  as  in  RM  i.e.,  both  are  held  constant  throughout 
the  simulation. 

Background  wind  and  temperature  fields  are  generated  prior  to 
starting  the  storm  simulation.  The  global  mean  heating  rate  Qq  is  just 
large  enough  to  maintain  the  initial  temperature  field  (Figure  2.1) 
from  cooling  off  due  to  thermal  conduction  (see  also  RM) .  In  addition, 
the  model  incorporates  EUV  and  UV  solar  heating  appropriate  to  the  sol¬ 
stitial  season  and  quiet-time  high-latitude  heating  appropriate  to 
solar  maximum.  Following  Dickinson  et.  al.  (1977),  the  magnitude  of 


12 


the  high-latitude  heating  was  adjusted  to  be  2.5  times  as  large  in  the 
summer  hemisphere  as  in  the  winter  hemisphere.  Figure  2.2  shows  the 
distribution  of  Q5  which  is  the  departure  from  the  global  mean  heating 
rate.  The  hemispherical  asymmetry  in  Qs  arises  from  the  greater  amount 
of  sunlight  in  the  summer  hemisphere  as  well  as  from  the  assumed  asym¬ 
metry  in  high-latitude  heating  just  mentioned. 

In  quantitative  terms,  Qs  is  approximately  80%  of  the  values  used 
by  Dickinson  et.  al.  (1977)  .  (See  Figure  9b  of  Dickinson  et .  al., 
p .  189 . )  It  should  be  noted  that  Qs  provides  the  thermal  drive  for  the 
quiet-time  solstitial  circulation. 

Using  Qs,  we  run  the  system  to  a  steady-state.  This  takes  approx¬ 
imately  5.5  days  of  simulation  time.  Figure  2.3  shows  the  time  histo¬ 
ries  of:  a)  the  temperature  change,  b)  the  meridional  winds,  and  c) 
the  zonal  winds  at  three  altitudes  at  40*  colatitude,  as  these  parame¬ 
ters  approach  a  steady-state.  It  can  be  seen  that  both  the  temperature 
and  the  meridional  winds  reach  steady-state  comparatively  quickly 
(within  2  to  3  days)  at  all  altitudes  as  opposed  to  the  zonal  winds 
which  take  between  5  to  6  days  at  the  lowest  altitudes  shown. 

Figure  2.4  shows  the  background  temperatures  and  wind  fields  at 
steady-state.  Figure  2.4a  shows  the  steady-state  temperature  field 
(solid  contour)  and  the  difference  temperature  field  (dashed  contour), 
which  is  the  difference  between  the  global  mean  temperature  (Figure 
2.1)  and  the  steady-state  temperature  field.  The  steady-state  tempe¬ 
ratures  are  a  result  of  a  balance  between  solar  and  high-latitude  heat¬ 
ing,  adiabatic  warming,  and  thermal  conduction.  The  difference  tempe- 
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rature  field  (dashed  contour)  shows  warming  by  as  much  as  200 ‘K  above 
150  km  in  the  summer  hemisphere.  However,  below  150  km,  there  is  a 
dipole-like  pattern  of  temperature  change  with  a  maximum  decrease  of 
40*  K  between  150  km  and  110  km,  and  a  maximum  increase  of  40*  K  be¬ 
low.  A  dipole-like  pattern  of  the  opposite  sense  exists  in  the  winter 
hemisphere  with  a  change  of  ±60*  K. 

Figure  2.4b  shows  a  streamline  plot  of  the  meridional  wind  field. 
The  plot  shows  rising  air  at  high  latitudes  in  the  summer  hemisphere. 
The  flow  proceeds  nearly  horizontal  towards  the  winter  hemisphere, 
crossing  the  equator,  and  sinking  at  approximately  140* -150*  colati¬ 
tude.  The  maximum  wind  speed  of  130  m/s  is  reached  at  30*  colatitude 
at  the  highest  altitudes.  A  reversed  circulation  also  exists  at  high 
latitudes  and  altitudes  in  the  winter  hemisphere.  It  extends  from  the 
pole  (180*)  to  140-150*  colatitude  and  exhibits  a  maximum  wind  speed 
of  15  m/s.  Below  120  km,  the  winds  blow  directly  from  summer  to  win¬ 
ter  hemispheres. 

Figure  2.4c  shows  the  zonal  winds.  Below  150  km,  a  two-celled  pat¬ 
tern  develops  with  maximum  winds  of  60-70  m/s,  westward  in  the  summer 
hemisphere  and  eastward  in  the  winter  hemisphere.  Above  150  km,  the 
zonal  winds  are  considerably  weaker.  In  this  region  the  ion  drag  force 
dominates  the  Coriolis  force  in  the  wind  balance.  The  summer  hemi¬ 
sphere  winds  above  150  km  maximize  at  30  m/s  in  the  westward  direction 

Our  winds  for  the  quiet-time  solstitial  circulation  are  similar  to 
the  earlier  two-dimensional  model  results  of  Dickinson  et.  al.  (1977) 
as  well  as  the  later  zonally  averaged  results  from  the  three-dimen- 
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sior.al  model  results  of  Dickinson  et .  al.  (1980).  However,  the  magni¬ 
tudes  of  our  winds  are  smaller,  owing  to  the  fact  that  our  values  of 
are  also  somewhat  smaller. 

"  STORM  SIMULATION  CONDITIONS 

As  in  Blanc  and  Richmond  (1980),  it  is  assumed  that  the  major 
f  ;urce  of  heat  input  for  our  thermospheric  storm  was  due  to  the  Joule 
dissipation  of  auroral  electro jet  currents.  Joule  dissipation  occurs 
when  a  current  flows  perpendicular  to  the  geomagnetic  field,  but  par¬ 
allel  to  the  electric  field.  Such  currents  are  known  as  Pedersen  cur¬ 
rents  and  the  corresponding  component  of  the  ionospheric  electrical 
conductivity  is  known  as  the  Pedersen  conductivity,  <7^  (Rishbeth  and 
Garriot,  1969) . 

In  deriving  a  reasonable  storm-time  heat  input  model  for  our  simu¬ 
lations,  we  have  consulted  two  geomagnetic  storms,  occurring  on  22 
March  1979  and  13  March  1989,  respectively.  These  storms  were  chosen 
because  they  showed  evidence  of  oscillations  in  the  auroral  activity 
index  (AE)  for  an  extended  period.  The  energy  inputs  for  these  storms 
(Figure  2.5a)  were  obtained  from  an  empirical  relation  (  Mazaudier  et . 
al.,  1987)  between  the  hemispheric  Joule  heating  and  the  auroral  acti¬ 
vity  index  (AE) .  In  this  manner,  it  was  found  'that  the  energy  inputs 

.  10  12 

varied  between  10  and  10  Watts  with  periods  of  2-4  hours  (typical 
periods  of  substorms) .  The  22  March  1979  storm  lasted  for  8  hours  and 
the  13  March  1989  storm  for  36  hours. 

Based  on  the  characteristics  of  the  two  storms  mentioned  above, 
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the  model  storm  energy  i^u  was  chosen  as  shown  in  Figure  2.5b.  The 

model  storm  lasts  for  12  hours.  It  shows  four  oscillations  with  peri- 

1 1 

cds  of  2-4  hours  and  peak  values  of  8.0  X  10  Watts  xn  energy  input. 

Our  model  for  the  colat itudinal  distribution  of  storm-energy  in¬ 
puts  is  also  based  on  available  information  on  Joule  dissipation  over 
the  polar  regions  for  typical  storms  (Ahn  et.  al.,  1983;  Kamide  et . 
al.,  1986;  Knipp,  1989).  Figure  2.6a  is  an  example  of  this  distribu¬ 
tion  (from  Knipp,  1989),  showing  that  most  of  this  dissipation  occurs 
in  a  narrow  region  of  colatitudes  centered  on  the  auroral  electrojet 
currents.  We  have  incorporated  this  feature  in  our  model  by  selecting 
the  colatitudinal  heating  profiles  to  be  Gaussian  in  shape,  having 
widths  of  10*  in  colatitude  and  centered  at  21*  and  159*  colatitude, 
repectivey  (Figure  2.6b). 

The  altitude  profile  for  the  storm-energy  input  (Figure  2.7)  is 
essentially  that  of  CJ^/p  (where  is  the  Pedersen  conductivity  )  which 
when  multiplied  by  ,  the  square  of  the  electric  field  assumed  to  be 
constant  with  height,  gives  the  heating  rate  per  unit  mass.  Because  of 
the  higher  air  density  at  lower  altitudes,  it  is  found  that  more  than 
80%  of  the  thermal  drive  for  the  thermospheric  storm-time  circulation 
is  provided  by  the  lower  maximum  in  this  profile.  Notice  that,  for 
reasons  of  simplicity  the  high-latitude  heat  inputs  for  the  storm  are 
assumed  to  be  hemispherically  symmetric. 
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FIGURE  CAPTIONS  FOR  CHAPTER  II 


Fig.  2.1.  The  global  mean  temperature  Tq  (*K). 

Fig.  2.2.  Qg  (10“®  Watts/Kg)  is  the  departure  from  the  global  mean 
heating  rate  of  the  solar  heating  at  EUV  and  UV  wavelengths  and  quiet¬ 
time  high-latitude  heating  at  solar  maximum. 

Fig.  2.3.  Time  histories  of:  a)  Temperature  changes  (*K),  b)  Meri¬ 
dional  winds  (m/s)  ,  and  c)  Zonal  winds  (m/s)  as  the  system  is  run  to 
a  steady-state.  These  time  histories  are  shown  at  altitudes  of  250  km 
(top  panel),  150  km  (middle  panel),  and  100  km  (bottom  panel)  at  40* 
colatitude. 

Fig.  2.4.  Steady-state  solstitial  background  temperature  and  wind 
fields : 

a)  Steady-state  temperature  (solid  line)  and  its  difference  from  the 
global  mean  temperature  (dashed  line)(*K),  b)  Background  meridional 
winds  with  streamlines  (m/s) ,  c)  Background  zonal  winds  (m/s)  with  E 
for  eastward  and  W  for  westward.  For  the  meridional  winds  the  barb 
notation  is  as  follows:  short  barb  (5  m/s);  long  barb  (10  m/s);  and 
flag  (50  m/s).  This  notation  is  also  followed  in  subsequent  figures. 

Fig.  2.5.  a)  The  estimated  hemispheric  energy  inputs  (in  units  of 
1011  Watts)  for  the  geomagnetic  storms  of  22  March  1979  (solid  line) 
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and  13-14  March  1989.  For  the  second  storm  lasting  for  36  hours,  the 
inputs  for  the  first  12  hours  (dashed  line)  and  the  second  12  hours 
(dotted  line)  are  shown  separately.  b)  Hemispheric  energy  input  for 
the  model  storm  (in  units  of  1011  Watts) . 

Fig.  2.6.  a)  Example  of  the  colatitudinal  distribution  of  Joule  heat¬ 
ing  from  Knipp  (1989),  and  b)  Normalized  colatitude  heating  profile  for 
the  model  strom. 

Fig.  2.7.  Normalized  altitude  heating  profile  (per  unit  mass)  for 
the  model  storm. 
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CHAPTER  III 


STORM  SIMULATION  RESULTS 

3.1  EVOLUTION  OF  STORM-TIME  WINDS  AND  TEMPERATURES 

The  onset  of  a  geomagnetic  storm  results  in  the  sudden  enhancement 
of  auroral  currents.  The  Joule  dissipation  of  these  currents  leads  to 
both  impulsive  heating  which  generates  equatorward  propagating  gravity 
wave  disturbances  in  both  summer  and  winter  hemispheres  as  well  as 
longer-lasting  heating  that  drives  a  global  meridional  circulation. 
The  global  circulation  arises  from  the  raising  of  the  geopotential 
heights  in  the  source  regions  resulting  in  a  pole-to-equator  pressure 
gradient  in  each  hemisphere.  In  this  chapter,  we  will  describe  these 
two  features  of  the  model  response  to  a  geomagnetic  storm. 

Figure  3.1  shows  the  storm-time  winds  at  2,2,1,  and  11  hours  after 
storm  onset.  At  2  hours  (Figure  3.1a)  the  earliest  wind  disturbances 
have  travelled,  at  high  altitudes,  a  distance  of  approximately  30*  in 
colatitude  from  the  regions  of  high-latitude  heating  in  either  hemi¬ 
sphere.  Above  250  km,  in  the  region  of  high  viscosity,  the  phases  of 
the  disturbances  are  uniform  with  height.  Below  250  km,  the  phase 
tilt  with  height  is  clearly  seen,  the  direction  of  the  tilt  being  per¬ 
pendicular  to  the  direction  of  wave  energy  propagation  which  is  upward 
and  away  from  the  source  region  (Hines,  1960;  Holton,  1975) .  In  the 
summer  hemisphere,  wind  speeds  associated  with  the  disturbances  exceed 
150  m/s  in  places,  which  is  approximately  50-100  m/s  above  the  back¬ 
ground  wind  speed  (Figure  2.3b) .  In  the  winter  hemisphere,  the  distur- 
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bance  wind  speeds  are  weaker. 

At  3  hours  (Figure  3.1b),  the  disturbances  from  both  hemispheres 
have  reached  the  equator.  At  7  hours  (Figure  3.1c)  the  disturbances 
from  the  summer  hemisphere  have  penetrated  well  into  the  winter  hemi¬ 
sphere.  Large-scale  meridional  circulation  patterns  are  also  evident, 
with  the  summer  hemisphere  cell  extending  to  140*  colatitude.  This 
cell  is  greatly  perturbed  by  the  disturbances  above  150  km.  The  winter 
hemisphere  cell  is  reversed  in  direction,  extending  to  110*  colatitude 
at  120  km  and  to  140*  colatitude  at  higher  altitudes. 

At  11  hours  (Figure  3 . Id) ,  the  disturbances  still  persist  above 
200  km,  but  below  this  altitude  the  meridional  circulations  have  be¬ 
come  better  defined.  The  summer  hemisphere  cell  meets  the  winter  hemi¬ 
sphere  ceil  at  approximately  110*  colatitude.  The  circulation  patterns 
continue  to  change  in  response  to  the  rapidly  varying  heat  sources 
assumed  in  the  model. 

Figure  3.2  features  the  storm-time  temperature  field  at  times  cor¬ 
responding  to  the  wind  fields  in  Figure  3.1.  At  2  hours,  above  200  km, 
temperatures  have  generally  increased  at  high-latitudes,  in  direct  re¬ 
sponse  to  storm-time  heating  (Figure  3.2a).  More  importantly,  the  fi¬ 
gure  shows  temperature  perturbations  of  50 *K  at  40-50*  colatitude  and 
120-140*  colatitude,  respectively.  Comparing  Figures  3.1a  and  3.2a, 
it  can  be  seen  that  the  temperature  perturbations  are  approximately  90* 
out  of  phase  with  the  vertical  velocities.  In  the  equatorial  region 
between  60*  and  140*  colatitudes,  the  temperature  field  is  undisturbed 
at  this  time. 
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At  3  hours,  however,  the  temperature  oscillations  have  reached 
the  equator  (Figure  3.2b),  and  mirror  the  vertical  velocity  oscilla¬ 
tions  (Figure  3.1b)  with  the  aforementioned  phase  shift.  At  the  same 
time,  at  high  latitudes,  evidence  for  adiabatic  cooling  resulting  from 
strong  upward  motion  above  the  heat  source  can  be  seen  in  the  decrease 
of  temperature  above  200  km. 

At  7  hours  (Figure  3.2c)  the  temperature  perturbations  have  cross¬ 
ed  over  the  equator  as  was  the  case  for  the  wind  perturbations.  Below 
250  km,  there  is  an  increase  in  temperature  on  the  order  of  100 *K  at 
mid  and  low  latitudes.  This  increase  is  apparently  associated  with 
the  establishment  of  the  meridional  circulation  (Figure  3.1c). 

At  11  hours  (Figure  3. 2d)  we  see  similar  features  in  the  tempera¬ 
ture  field  as  at  7  hours  below  250  km.  Above  this  height,  additional 
regions  of  warming  of  50-100 *K  due  to  the  persistence  of  wind  distur¬ 
bances  are  also  observed. 

3.2  INFLUENCE  OF  BACKGROUND  SEASONAL  WINDS  AND  TEMPERATURES  ON  MEAN 

STORM-TIME  WINDS  AND  TEMPERATURES 

Figure  3.3  shows  the  12-hour-averaged  storm-time  wind  fields  for 
model  runs  with  and  without  background  seasonal  winds  and  associated 
temperatures  (Figures  3.3a  and  b) .  The  averaged  wind  fields  are  ob¬ 
tained  from  the  instantaneous  wind  fields  at  10  minute  intervals  over 
the  12  hour  span  of  our  model  storm. 

With  background  winds,  the  storm- time  meridional  circulation  shows 
a  two-celled  structure  with  a  larger  summer  hemisphere  cell  extending 
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to  110-130*  colatitude.  Both  cells  show  increases  in  wind  speed  over 
the  background  solstice  winds  (Figure  2.3b).  This  is  best  shown  by 
the  reversed  winter  hemisphere  cell  which  exhibits  winds  averaging  50 
m/s  compared  to  10-15  m/s  for  the  background  winds  alone.  The  latitu¬ 
dinal  extent  of  this  cell  is  also  greater,  extending  to  110*  colatitude 
at  150  km.  In  addition,  the  vertical  extent  of  this  cell  has  increased 
and  the  center  of  the  cell  has  shifted  downward  from  175  km  to  120  km. 

In  comparison,  the  storm  circulation  without  background  winds  shows 
a  two-celled  structure,  with  the  two  cells  being  symmetric  about  the 
equator.  This  symmetry  merely  mirrors  the  corresponding  symmetry  in 
the  heat  sources  (Figure  2.6b). 

The  zonal  winds  generated  by  the  storm  consist  of  strong  westward 
jets  equatorward  of  the  heat  sources  and  much  weaker  eastward  flows 
poleward  (Figure  3.3d).  This  storm-time  zonal  wind  system  enhances 
the  background  winds,  as  can  be  seen  by  comparison  of  Figures  3.3c  and 
2.3c.  This  enhancement  leads  to  the  intensification  and  poleward  mi¬ 
gration  of  the  westward  jet  normally  present  in  the  summer  hemisphere 
during  the  solstitial  season.  In  the  winter  hemisphere,  we  have  a 
regime  with  a  storm-time  westward  jet  located  overhead  of  the  seasonal 
eastward  jet  with  strong  shears  in  the  intermediate  zone. 

Lacking  direct  observational  support,  it  is  not  clear  to  what 
extent  the  double-jet  structure  of  the  zonal  wind  system  and  the  inter¬ 
mediate  region  of  strong  shears  calculated  for  the  winter  hemisphere 
are,  in  fact,  realistic.  It  should  also  be  recalled  that  there  are 
additional  storm-time  winds  generated  by  ion  drag  effects  near  the 
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auroral  electrojet  region  which  are  not  included  in  our  model 
(Hernandez  and  Roble,  1976;  Killeen  et .  al.  1984,  Killeen  et .  al., 
1986)  . 

Figure  3.4  shows  corresponding  results  for  the  12-hour-averaged 
storm-time  temperature  fields  with  background  winds  (Figure  3.4a)  and 
without  background  winds  (Figure  3.4b).  Careful  inspection  shows  that 
for  both  these  cases,  there  is  an  average  increase  of  50*  K  in  the 
storm-time  temperatures  at  middle  and  low  latitudes.  This  temperature 
increase  is  due  to  downward  velocities  of  the  storm-generated  part  of 
the  meridional  circulation,  whose  strength  is  nearly  the  same  with  or 
without  background  winds. 


3.3  CHARACTERISTICS  OF  STORM-GENERATED  WAVE  DISTURBANCES 

We  next  turn  our  attention  to  the  disturbances  generated  by  the 
storm.  Figures  3.5a  and  3.5b  exhibit  the  storm-generated  disturbances 
at  two  colatitudes  at  an  altitude  of  325  km.  The  top  panel  of  Figure 
3.5a  shows  the  heat  input  profile  (in  units  of  101 1  Watts),  the  middle 
panel  shows  the  vertical  velocity  component  of  the  wind  response  at 
40*  colatitude,  and  the  bottom  panel  shows  the  vertical  velocity  re¬ 
sponse  at  70*  colatitude.  In  Figure  3.5b,  the  meridional  component 
of  the  wind  response  is  shown  instead  of  the  vertical  velocity  compo¬ 
nent  . 

The  oscillating  heat  source  generates  wave  disturbances  throughout 
the  12-hour  period  of  storm  simulation.  We  define  a  disturbance,  some¬ 
what  arbitrarily,  as  a  full  cycle  of  the  oscillation  in  vertical  velo- 
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city.  Six  disturbances  are  identified  in  this  manner  and  are  numbered 
in  Figures  3.5a  and  b.  The  beginning  and  ending  times  of  each  distur¬ 
bance  are  also  indicated  by  vertical  markers.  At  40*  colatitude  these 
times  are  0:30-2:00  storm-time  (ST)  for  the  first  disturbance,  2:00- 
4:00  ST  for  the  second,  4:00-6:00  ST  for  the  third,  6:00-7:50  ST  for 
the  fourth,  7:50-10:00  ST  for  the  fifth,  and  finally  10:00-12:00  ST 
for  the  sixth  disturbance. 

The  first  disturbance,  which  arrives  at  40*  colatitude  as  early  as 
0:30  ST,  must  have  originated  from  the  first  heat  pulse  in  the  form  of 
a  spectrum  of  gravity  waves  with  speeds  up  to  900  m/s.  The  initial 
part  of  the  second  disturbance  before  2:30  ST  must  also  have  originat¬ 
ed  from  the  comparatively  slower  gravity  waves  generated  by  the  first 
heat  pulse.  At  2:30  ST  the  fastest  wave  from  the  second  heat  pulse 
should  have  arrived  at  this  colatitude,  if  its  velocity  was  the  same 
(-900  m/s)  as  of  the  first  wave  from  the  initial  pulse.  Hence,  waves 
from  both  the  first  and  second  heat  pulses  would  contribute  to  the 
subsequent  portions  of  the  second  disturbance.  The  most  probable  in¬ 
terpretation  for  the  third  disturbance  between  4:00  and  6:00  ST  is 
that  it  results  from  the  slower  waves  from  both  the  first  and  second 
heat  pulses.  Similar  to  the  first  two  disturbances,  the  fourth  and 
fifth  disturbances  consist  of  waves  mostly  from  the  third  and  fourth 
heat  pulses,  respectively,  and  to  a  lesser  extent,  from  the  earlier 
pulses.  The  interpretation  of  the  sixth  disturbance  is  similar  in 
nature  to  that  for  the  third  disturbance. 

In  the  bottom. panel  of  Figure  3.5a  for  70*  colatitude,  much  of  the 
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finer  structure  in  the  disturbances  has  disappeared.  Comparison  of 
the  middle  and  bottom  panels  in  this  figure  show  that  the  disturbances 
travel  from  40*  to  70*  colatitude  with  a  time  delay  of  70-120  minutes. 

The  meridional  velocity  disturbances  in  Figure  3.5b,  whose 
durations  are  assumed  to  be  the  same  as  the  corresponding  vertical 
velocity  disturbances  in  Figure  3.5a,  have  certain  peculiar  features. 
For  instance,  the  latter  part  of  the  second  disturbance  after  5:00  ST 
and  all  of  the  third  disturbance  at  70*  colatitude  are  considerably 
depressed  in  magnitude  from  their  magnitudes  at  40*  colatitude.  In 
addition,  the  fourth  disturbance  (6:00-7:50  ST)  and  the  earlier  part 
of  the  fifth  disturbance  (8:00-8:50  ST)  passing  40*  colatitude  are 
depressed  in  magnitude  compared  to  their  magnitudes  at  70*  colatitude'. 
These  features  are  most  likely  explained  as  due  to  a  disturbance 
launched  in  the  southern  hemisphere,  which  is  propagating  in  the  oppo¬ 
site  direction  to  the  other  disturbances.  As  this  southern  hemispheric 
disturbance  propagates  from  70*  to  40*  colatitude  with  a  speed  of  900 
m/s,  it  is  successively  superposed  on  the  equatorward  propagating  dis¬ 
turbances  causing  the  ’’depression”  in  magnitude  of  the  second  and 
third  disturbances  at  70*  colatitude  and  later  the  fourth  and  fifth 
disturbances  at  40*  colatitude. 

Figure  3.5a  also  gives  information  concerning  the  phase  speeds, 
dispersion  and  dissipation  characteristics  of  the  disturbances.  The 
phase  speed  of  the  leading  edge  of  the  first  disturbance  is  approxi¬ 
mately  790  m/s  as  it  propagates  from  40*  to  70*  colatitude.  Other 
wave  disturbances  in  Figure  3.5a  exhibit  slower  phase  speeds  averaging 
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450-550  m/s ,  which  is  in  line  with  the  observations  of  Hernandez  and 
Roble  (1978)  for  storm-generated  disturbances  at  middle  latitudes. 
We  also  see  evidence  of  dispersion  by  the  increase  of  the  periods  of 
the  disturbances  as  they  travel  from  40*  to  70*  colatitude.  For  exam¬ 
ple,  the  first  disturbance  at  40*  colatitude  has  a  period  of  90  min¬ 
utes.  When  this  disturbance  reaches  70*  colatitude,  its  period  has 
increased  to  120  minutes.  Similarly,  the  second  disturbance  has  a 
period  of  120  minutes  at  40*  colatitude  and  a  period  of  145  minutes  at 
70*  colatitude.  Furthermore,  Figure  3.5a  also  shows  that  fast-moving, 
long-period  waves  tend  to  maintain  their  identities  over  longer  dis¬ 
tances,  as  compared  to  slower-moving  high-frequency  waves  because  they 
are  less  affected  by  dissipation  (Hines,  1960;  Richmond,  1978a)  .  This 
can  be  seen  most  clearly  in  the  longer  period  of  the  disturbances  in 
regions  further  away  from  the  source  region. 

3.4  INFLUENCE  OF  BACKGROUND  WINDS  AND  TEMPERATURES 

The  influence  of  the  background  meridional  winds  and  temperatures 
on  the  wave  disturbances  is  illustrated  in  Figure  3.6,  where  we  observe 
velocity  disturbances  that  started  in  the  summer  and  winter  hemispheric 
source  regions  and  reached  40*  (solid  line)  and  140*  (dashed  line)  co¬ 
latitudes,  respectively,  having  travelled  the  same  meridional  distance 
In  the  summer  hemisphere,  the  disturbances  are  propagating  in 
the  direction  of  the  background  meridional  wind  (Figure  3.7).  In  the 
winter  hemisphere,  the  disturbance  propagation  is  opposite  to  the  di¬ 
rection  of  meridional  winds.  In  addition,  the  disturbances  in  the 
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summer  hemisphere  are  propagating  in  a  region  of  warmer  temperatures 
than  the  disturbances  in  the  winter  hemisphere.  As  a  result,  the  dis¬ 
turbances  in  the  summer  hemisphere  are  observed  to  arrive  at  40*  cola¬ 
titude,  earlier  than  the  arrival  of  the  disturbances  in  the  winter 
hemisphere  at  140*  colatitude.  This  time  difference,  however,  amounts 
to  no  more  than  about  12  minutes  in  a  travel  time  of  nearly  an  hour, 
even  at  the  highest  altitude  (325  km)  shown  in  Figure  3.6. 

In  order  to  understand  the  observed  time  difference,  mentioned 
above,  one  may  start  with  the  simplified  expression  for  the  intrinsic 
frequency  for  a  gravity  wave  in  a  hydrostatic  and  isothermal  atmo¬ 
sphere.  For  our  log  p  coordinate  system  (Chapter  2)  with  Z=ln(ps/p), 
this  frequency  is  given  by: 

W2  =  N2  H2  (  ~-1  )  (3.1) 

“  +  4 

where  N  is  the  buoyancy  frequency,  H  is  the  scale-height,  and  Z  and  m 
are  the  meridional  and  vertical  wave  numbers,  respectively.  The  above 
expression  for  the  intrinsic  frequency  is  the  same  as  that  found  in 
Andrews  et .  al.  (1987)  provided  we  note  that  in  the  log  p  coordinate 

system  used  by  these  authors  Z=H  In(ps/p) . 

2  2 

Using  the  isothermal  value  N  =  Kg/H  =  KRT/H  ,  where  K=R/Cpf  we 
obtain  the  intrinsic  phase  velocity  in  the  meridional  direction  as: 
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which  shows  that  the  intrinsic  phase  velocity  varies  directly  as  the 
square  root  of  the  temperature. 

We  observe  from  Figure  2.4a  that  the  average  temperature  at  325  km 
between  the  source  region  and  40*  colatitude  in  the  summer  hemisphere 
is  1080*  K.  Similarly,  the  average  temperature  between  the  source 
region  and  140*  colatitude  in  the  winter  hemisphere  is  960*  K.  Using 
these  temperatures  values  in  the  above  expression  for  the  phase  velo¬ 
city,  we  find  that  the  ratio  of  the  phase  velocities  of  the  waves 
travelling  equal  distances  in  the  two  hemispheres  is  in  the  ratio 
(1080/960) =  1.06,  accounting  for  a  3.6  minute  difference  in  travel 
times  in  a  period  of  one  hour  due  to  the  temperature  difference 
between  the  hemispheres. 

In  addition,  the  intrinsic  phase  velocities  of  the  waves  (-500  m/s) 
are  Doppler-shif ted  by  the  background  winds.  This  Doppler-shif t  is 
positive  in  the  summer  hemisphere  (where  the  background  meridional 
winds,  as  pointed  out  earlier,  are  in  the  same  direction  as  the  direc¬ 
tion  of  wave  propagation  (equatorward) ,  and  have  an  average  speed  of 
100  m/s  along  the  wave  path  at  325  km.  In  the  winter  hemisphere,  the 
wind  direction  changes  along  the  wave  path,  avey.v  ing  out  to  almost 
zero.  Because  of  the  Doppler-shift  effect,  the  resultant  phase  velo¬ 
cities  of  the  waves  in  the  summer  and  winter  hemispheres  are,  there¬ 
fore,  in  the  ratio  600/500  =  1.2,  accounting  for  a  12  minute  differ¬ 
ence  in  a  travel  time  of  one  hour. 

Our  simple  calculation,  therefore,  gives  a  total  sixteen  minute 
time  difference  in  the  travel  times  of  the  disturbances  at  325  km  as 
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against  the  twelve  minute  difference  more  accurately  calculated  in  our 
model.  At  lower  altitudes,  where  most  of  the  wave  energy  resides, 
both  the  wind  and  temperature  effects  on  phase  velocities  and  travel 
times  become  negligible. 

Figure  3.6  also  shows  that  the  amplitudes  of  the  disturbances  in 
the  summer  hemisphere  travelling  in  the  direction  of  the  background 
meridional  wind  are  less  than  for  disturbances  in  the  winter  hemisphere 
travelling  against  the  background  flow.  This  is  to  be  expected  based 
on  the  approximate  tendency  for  the  waves  in  the  disturbances  to  con¬ 
serve  wave  action  density  A=E/  (CO-£V)  ,  where  E  is  the  wave  energy  den¬ 
sity,  (0  is  the  intrinsic  frequency,  i  is  the  wave  number,  V  in  the 
wind  in  the  direction  of  wave  propagation,  and  £- kV  is  the  Doppler- 
shifted  frequency  (Andrews  et.  al.,  1987;  Bretherton  and  Garrett, 
1968)  .  As  0 >-£V  decreases,  as  in  the  example  of  summer  hemispheric 
disturbances,  the  wave  energy  (and  hence  amplitude)  must  decrease,  as 
is  observed.  The  net  effect  is  small  because  the  phase  speeds  of  the 
disturbances  are  much  larger  then  the  background  wind  velocity. 

3.5  VERTICAL  STRUCTURE  OF  THE  WAVE  DISTURBANCES 

As  illustrative  of  the  vertical  structure  of  the  disturbances, 
Figure  3.8a  shows  the  first  two  vertical  velocity  disturbances  (solid 
line)  and  the  corresponding  temperature  disturbances  (dashed  line)  at 
70*  colatitude  for  altitudes  between  100  and  450  km.  For  the  velocity 
disturbances,  the  distance  between  two  consecutive  altitude  markers 
also  represents  a  vertical  velocity  of  5  m/s.  The  maximum  amplitude 
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of  the  velocity  disturbances  is  ±7.5m/s.  For  the  temperature  distur¬ 
bances  the  same  distance  between  the  altitude  markers  corresponds  to 
50  *  K .  The  maximum  amplitude  of  the  temperature  disturbances  is  50*K. 
The  amplitudes  of  the  velocity  disturbances  do  not  increase  as  rapidly 
as  the  inverse-root  of  air  density  because  of  the  increasing  impor¬ 
tance  of  viscous  dissipation  with  height.  Increasing  dissipation  is 
also  responsible  for  the  decrease  of  phase  tilt  with  height  (Hines, 
1968),  the  phase  becoming  uniform  above  250  km  (see  also  Figure  3.1). 
We  further  notice  that  the  phase  difference  between  the  vertical  velo¬ 
city  disturbances  and  the  temperature  disturbances  is  close  to  90*, 
as  would  be  the  case  if  adiabatic  heating  (cooling)  due  to  the  verti¬ 
cal  motions  of  the  velocity  disturbances  is  the  cause  of  the  tempera¬ 
ture  disturbances . 

Figure  3.8b  shows  the  same  disturbances  in  meridional  velocity 
(solid  line)  with  the  distance  between  the  altitude  markers  equal  to 
50  m/s  (positive  equatorward) .  The  meridional  velocity  disturbances 
show  a  similar  increase  of  amplitude  with  height  and  a  similar  phase 
tilt  as  the  vertical  velocity  disturbances.  The  meridional  velocity 
disturbances  are  different,  however,  in  that  they  are  nearly  in  phase 
with  the  temperature  disturbances.  This  means  chat  the  vertical  and 
meridional  velocities  are  out  of  phase  by  90*,  which  follows  from  the 
well-know  polarization  relations  (Gossard  and  Hooke,  1977/  Hines, 
1960)  . 

3.6  WAVE-INDUCED  HEAT  AND  MOMENTUM  FLUXES 
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Both  the  vertical  and  meridional  velocity  disturbances  can,  over 
the  period  of  the  disturbances,  produce  net  fluxes  of  heat  and  meri¬ 
dional  momentum  in  the  vertical  and  meridional  directions,  respective¬ 
ly.  For  example,  a  net  heat  flux  is  produced  if  the  vertical  (meridi¬ 
onal)  velocity  disturbances  and  the  temperature  disturbances  are  not 
in  exact  quadrature  (i.e.  phase  difference  is  not  90*)  or  if  the  dis- 
bances  are  not  exactly  symmetric  (i.e.,  the  average  velocity  or  the 
net  temperature  change  over  the  period  of  the  disturbance  is  not  zero) 
either  due  to  dissipation  or  non-stationarity  (i.e.,  change  in  wave 
amplitude  over  a  wave  period) . 

Figures  3.9a,b,  and  c  show  the  vertical  and  meridional  heat  fluxes, 
and  the  vertical  flux  of  meridional  momentum  generated  at  250  km  by 
the  disturbances  shown  in  Figures  3.8a  and  3.8b,  Figure  3.9a  shows 
the  vertical  velocity  disturbances  (top) ,  the  temperature  disturbances 
(middle) ,  and  the  vertical  heat  fluxes  (bottom) .  The  starting  and 
ending  times  of  the  disturbances  to  be  referred  to  below,  are  with 
respect  to  the  vertical  velocity  disturbance  as  observed  at  the  spec¬ 
ific  altitude  of  250  km.  The  first  vertical  velocity  wave  (1:40-4:00 
ST)  is  nearly  in  quadrature  with  its  corresponding  temperature  wave. 
In  addition,  it  is  only  slightly  asymmetric  in  shape.  Hence,  its  net 
vertical  heat  flux  appears  to  be  close  to  zero.  The  second  vertical 
velocity  wave  (4:00-6:15  ST)  is  also  nearly  in  quadrature  with  its  cor¬ 
responding  temperature  wave.  However,  this  wave  is  more  asymmetric 
in  shape  (i.e.  the  downward  velocities  are  larger  then  the  upward  velo¬ 
cities)  and  yields  a  net  downward  flux  of  heat. 
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Correspondingly,  Figure  3.9b  shows  the  meridional  velocity  distur¬ 
bances  (top) ,  the  temperature  disturbances  (middle) ,  and  the  meridional 
heat  fluxes  (bottom).  As  previously  mentioned  (Chapter  3.3,  Figure 
3.5b)  during  the  period  (4:00-6:15  ST)  of  the  second  disturbance, 
there  is  indication  of  superposition  effects  from  a  poleward  travel¬ 
ling  disturbance  as  well.  Wha;  we  notice  is  that,  the  meridional  velo¬ 
city  disturbances  are  asymmetric  in  shape  although  in  phase  with  the 
corresponding  temperature  disturbances,  resulting  in  a  net  equatorward 
heat  flux. 

Similarly,  Figure  3.9c  shows  the  vertical  velocity  (top),  the 
meridional  velocity  (middle)  disturbances,  and  the  vertical  flux  of 
meridional  momentum  (bottom) .  The  first  vertical  velocity  disturbance 
and  the  corresponding  meridional  velocity  disturbance  (1:40-4:00  ST) 
are  neither  in  quadrature  with  each  other  nor  are  symmetric  in  shape. 
This  results  in  a  net  upward  meridional  momentum  flux. 

Between  4:00  and  6:15  ST,  when  the  second  disturbance  is  also 
influenced  by  a  poleward  travelling  disturbance  toward  the  end  of  its 
period,  it  is  the  downward  velocity  peak  that  is  closer  in  time  to  the 
positive  peak  in  the  meridional  velocity  and,  as  a  result,  the  net  ver¬ 
tical  flux  of  meridional  momentum  is  downward. 
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FIGURE  CAPTIONS  FOR  CHAPTER  III 


Fig.  3.1.  Time  evolution  of  storm-time  winds.  The  winds  (m/s)  are 
shown  at:  a)  2  hours,  b)  3  hours,  c)  7  hours,  and  d)  11  hours  after 
storm  onset. 

Fig.  3.2.  Same  as  Figure  3.1  except  for  storm-time  temperatures  (*K). 

Fig.  3.3.  Time-averaged  storm-time  meridional  winds:  a)  With  back¬ 
ground  winds  and  b)  Without  background  winds.  Time-averaged  storm¬ 
time  zonal  winds:  c)  With  background  winds  and  d)  Without  background 
winds.  (See  text  for  definitions  of  with  and  without  background 
winds . ) 

Fig.  3.4.  Same  as  Figure  3.3  except  for  storm-time  temperatures:  a) 
With  background  winds  and  b)  Without  background  winds. 

Fig.  3.5.  Storm-generated  wave  disturbances:  a)  Disturbances  in  ver¬ 
tical  velocity  at  325  km  are  shown  at  40*  colatitude  in  the  middle 
panel  and  at  70*  colatitude  in  the  bottom  panel.  The  top  panel  is 
the  storm-time  hemispheric  heat  input  (in  units  of  10^  Watts)  profile 
shown  for  the  purpose  of  relating  the  waves  to  their  sources,  b)  Same 
as  a)  except  for  disturbances  in  meridional  velocity. 

Fig.  3.6.  Influence  of  background  seasonal  winds  and  associated  tem- 
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peratures  on  vertical  velocity  disturbances  propagating  from  the  sum¬ 
mer  and  winter  hemispheric  source  regions  (see  text) .  The  disturbances 
at  40*  colatitude  (solid  line)  and  at  140*  colatitude  (dashed  line) 
are  shown  at  325  km  (top  panel)  ,  250  km  (middle  panel)  ,  and  150  km 
(bottom  panel) . 

Fig.  3.7.  Background  meridional  winds  for  the  solstitial  season.  The 
corresponding  zonal  winds  are  shown  in  Figure  2.4c. 

Fig.  3.8.  Vertical  structure  of  the  first  two  waves  at  70*  colatitude 
(see  text) :  a)  Vertical  velocity  (solid  line)  and  temperature  change 
(dashed  line) .  b)  Meridional  velocity  (solid  line)  and  temperature 
change  (dashed  line)  . 

Fig.  3.9.  Transport  of  heat  and  momentum  by  the  first  two  wave  distur¬ 
bances  at  250  km.  The  parameters  shown  in  this  figure  are:  a)  Verti¬ 
cal  velocity  (top  panel) ,  Temperature  change  (middle  panel) ,  and  Ver¬ 
tical  heat  flux  (bottom  panel).  b)  Meridional  velocity  (top  panel)/ 
Temperature  change  (middle  panel) /  and  Meridional  heat  flux  (bottom 
panel) .  c)  Vertical  velocity  (top  panel) f  Meridional  velocity  (middle 
panel) /  and  Vertical  flux  of  meridional  momentum  (bottom  panel) . 
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CHAPTER  IV 


WAVE-MEAN  FLOW  INTERACTION 


4.1  EFFECTS  ON  CIRCULATION 

To  investigate  the  effects  of  wave-mean  flow  interaction  for  our 
model  storm,  we  partition  the  dependent  variables  (V^  Vq^  w,  and  T) 
in  the  momentum  and  thermodynamic  energy  equations  into  a  time-averaged 
part  and  an  eddy  part,  defined  as  the  deviation  from  the  time-average. 
The  time-averaged  part  is  taken  to  be  the  storm-time  response,  averaged 
over  12  hours,  of  the  thermosphere  with  seasonal  background  winds  and 
temperatures.  (In  the  following  discussion,  we  use  the  words  eddies 
and  waves  interchangeably.)  This  partition  yields: 


II 
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v4> 
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<r> 

n 

<1 

CD 

+ 

v0 

w  =  w 

+ 

t 

w 

T  =  T  +  T  with  T  -  Tq  +  T^ 


(4.1) 


where  the  overbar  represents  the  time  average  and  the  prime  represents 
the  deviation  from  this  average.  In  addition,  is  the  difference 
between  the  time-averaged  (T)  and  initial  (Tq)  temperatures. 

We  then  substitute  the  dependent  variables  into  the  equations  of 
motion  (Eq.  2.1,  2.2,  and  2.4)  and  the  thermodynamic  energy  equation 
(Eq.  2.5),  and  average  over  a  12-hour  period.  With  slight  rearrange¬ 
ment  the  resulting  time-averaged  equations  are; 
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The  zonal  momentum  equation: 


dv6 

at 


ve  av*  _  av*  _  vavq 

TaT  -  waf  "  2%cos0  -  ~^cot0 


(4.2) 


is  h  ( )  -  i  fe  ( p-S )  - 


a  SJ 


1  3 
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v<t>  ve 


+  X 


The  meridional  momentum  equation: 
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i  a 


a  sin8  80 


a  r  '2  \  i  a  (  '  '  \  v<b 
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Continuity  equation: 


be  ae  (  vesin0  )  +  «  fe  (  pw  )  =  0 


a  sinG  30 


1  L 
p  3z 


(4.4) 


Thermodynamic  energy  equation: 


at  (cpT)  ~ 


ve  a 

"  a  ae  (cpt) 
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w  (CpT)  -  RwT 


(4.5) 


i  a 


a  sine  ae  (  cpveT  sin0  )  -  p  fi  cp(pw’t')  -  ^'T' 
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+  Q  +  K  . 

In  practice,  the  time-averaged  meridional  circulation  (Vq,  w  )  de¬ 
scribed  by  the  above  equations  was  obtained  as  the  average  at  fixed 
grid  points  of  the  instantaneous  solutions  of  the  model  equations. 
This  circulation,  also  referred  to  as  the  mean  Eulerian  circulation, 
is  shown  in  Figure  3.3a.  We  can  now  alternatively  consider  the  mean 
Eulerian  circulation  as  driven  by  the  forcing  terms  on  the  right-hand 
sides  of  the  equations  above.  In  doing  so,  it  is  convenient  to  separ¬ 
ate  the  linear  and  non-linear  forcing  terms,  and  think  of  these  terms 
as  driving  a  mean  linear  circulation  and  a  mean  induced  circulation, 
respectively . 

The  procedure  that  was  adopted  to  obtain  the  mean  linear  circula¬ 
tion  (Vq1,  w1)  was  to  run  the  model  with  storm-time  energy  inputs 
reduced  to  l/100th  of  their  original  strength.  It  was  assumed  that 
the  resulting  model  response  would  be  essentially  linear.  This  re¬ 
sponse  is  scaled  back  up  by  a  factor  of  100  and  is  shown  in  Figure 
4.1  as  the  mean  linear  circulation 

We  notice  that  the  mean  Eulerian  circulation  (Figure  3.3a)  and  the 
mean  linear  circulation  are  remarkably  similar  in  both  structure  and 
magnitude.  Both  circulations  exhibit  a  two-celled  structure  with  the 
larger  summer  hemisphere  cell  extending  across  the  equator,  meeting  a 
reversed  winter  cell  at  100-120*  colatitude,  and  resulting  in  a  region 
of  strong  downward  motion  there.  The  similarity  of  these  circulations 
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implies  that  the  mean  linear  circulation  accounts  for  a  major  part  of 
the  mean  Eulerian  circulation. 

The  difference  between  the  mean  Eulerian  and  mean  linear  circula- 

—  *  _ * 

tions  is  the  mean  induced  circulation  (Vq  ,  w  )  shown  in  Figure  4.2. 

The  mean  induced  circulation  results  from  the  non-linear  forcing  terms 

on  the  right-hand  sides  of  the  time-averaged  momentum  and  thermodynamic 

energy  equations  (4 . 2, 4 . 3/ 4 . 5) .  The  non-linear  terms  consist  of  two 

types:  those  involving  the  products  of  two  mean  quantities  (e.g. 

w  ^  Vq)  and  those  involving  the  products  of  two  eddy  quantities  (e.g. 
1  3  1 

—  0-^  (pw'Vg)  ) .  Detailed  calculations  show,  however,  that  the  eddy 
non-linear  terms  are  considerably  more  important  in  forcing  the  induced 
circulation  than  the  mean  non-linear  terms,  except  in  the  region  of 
storm-time  forcing.  Neglecting  the  mean  non-linear  terms,  the  govern¬ 
ing  equations  for  the  induced  circulation  are: 
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Figure  4.2  shews  that  the  most  prominent  features  of  the  induced 
circulation  are  the  two  cells  on  either  side  of  the  equator,  at  alti¬ 
tudes  below  200  km.  These  cells  have  the  same  sense  of  rotation  as 
the  storm-time  mean  Eulerian  circulation  cells  (Figure  3.3)  and  there¬ 
fore,  make  a  contribution  to  the  strengths  of  the  latter.  Above  200 
km,  the  induced  circulation  tends  to  split  into  three  or  more  cells 
of  more  limited  extent  in  colatitude.  It  is  important  to  remark,  in 
this  context,  that  the  above  mentioned  features  of  the  induced  circu¬ 
lation  do  not  substantially  change  if  we  adopt  somewhat  different 
averaging  periods  from  the  12  hour  period  chosen  here,  They  are, 
therefore,  robust  features  of  the  storm-time  dynamics. 

4.2  HEAT  SOURCES  AND  SINKS 

The  mean  Eulerian,  the  mean  linear,  and  the  mean  induced  circula¬ 
tions  set  up  their  own  heat  sources  and  sinks.  These  heat  sources  and 
sinks  generated  by  the  vertical  and  horizontal  motions  associated  with 
these  circulations  are  expressed,  in  what  follows,  as  equivalent  tem¬ 
perature  changes  and  are  compared  with  one  another. 

The  equivalent  temperature  change  (AT)  produced  by  the  mean 
Eulerian  circulation  (Vq,  w  )  over  a  period  At  =  12  hours  is: 

"fe  <cp?>  -  r  fe  (CP?I  )  <4-9> 

The  resulting  pattern  shown  in  Figure  4.3  exhibits  alternating  regions 
of  comparatively  large  and  small  temperature  changes,  forming  three 
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"couplets",  above  200  km  between  40*  and  140*  colatitudes.  For  exam¬ 
ple,  one  such  couplet  pattern  with  maximum  and  minimum  temperatures 
of  220*  K  and  160*  K,  repectively,  occur  over  the  equator. 

For  the  mean  linear  circulation,  the  equivalent  temperature  changes 
over  the  12-hour  period  can  be  similarly  expressed  as: 


AT 


At  ( 

cp  1 


-Rw^T  - 


1,1  k  icpt) 


A  <L 

a  30 


(4.10) 


The  resulting  pattern  in  Figure  4.4  shows  an  overall  correspondence  to 
temperature  changes  due  to  the  mean  Eulerian  circulation.  This  in¬ 
cludes  the  large  temperature  increase  associated  with  the  main  region 
of  descent  at  110-120*  colatitude.  However,  the  couplet  patterns  of 
temperature  changes  associated  with  the  mean  Eulerian  circulation  at 
higher  altitudes  are  not  reproduced  in  the  temperature  changes  due  to 
the  mean  linear  circulation. 

The  differences  between  Figures  4.3  and  4.4  are  the  equivalent 
temperature  changes  to  be  attributed  to  the  induced  circulation,  and 
is  shown  in  Figure  4.5.  This  difference  pattern  should  capture  the 
couplet  patterns  of  temperature  changes  associated  with  the  mean 
Eulerian  circulation.  The  couplet  patterns  are,  therefore,  judged  to 
be  the  signature  of  wave-mean  flow  interaction  on  the  temperature  field 
of  the  mean  Eulerian  circulation. 

The  equivalent  temperature  changes  due  to  the  induced  circulation 
can  now  be  compared  to  the  changes  produced  by  the  mean  eddy  heat  flux 
convergence  (divergence) : 
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AT 


At  r  '  '  1  d 

■  c;  l-R“T  • ;  Ji  cp<p"  T 1 


a  sin0  30 


30  (CpVQT’sinG))  (4.11) 


These  temperature  changes  are  shown  in  Figure  4.6. 

Comparison  of  Figures  4.5  and  4.6  shows  that,  at  least  in  middle 
latitudes,  regions  of  warming  (cooling)  due  to  the  induced  circulation 
roughly  correspond  to  regions  of  cooling  (warming)  due  to  mean  eddy 
heat  flux  divergence  (convergence) .  However,  this  correspondence  does 
not  hold  near  the  equator,  where  both  the  induced  circulation  and 
eddy  heat  flux  convergence  cooperate  to  produce  a  region  of  warming. 
This  equatorial  warming  is  on  the  order  of  50-90*  K  above  150  km,  and 
is  1/3  to  1/2  of  the  warming  produced  by  the  mean  linear  circulation. 
The  equatorial  warming  due  to  the  induced  circulation  alone  is  40- 
60*  K,  and  is  roughly  1/4  to  1/3  of  the  warming  due  to  the  mean 
Eulerian  circulation. 

It  is  important  to  emphasize  that  the  above  discussion  of  tempera¬ 
ture  changes  due  to  the  different  circulations  is  with  respect  to 
equivalent  and  not  actual  temperature  changes.  To  discuss  the  actual 
temperature  changes,  we  need  to  consider  the  dissipative  effects  of 
molecular  viscosity  and  thermal  conductivity.  Our  detailed  heat  bud¬ 
get  calculations  including  these  effects  are  not  given  here. 

The  non-cancellation  of  the  thermal  effects  of  the  induced  circu¬ 
lation  and  the  mean  eddy  heat  flux  convergence  (divergence)  is  due  to 
the  temporal  terms  (3VQ*/3t  ,  3T*/3t)  on  the  left-hand  side  of  the 
governing  equations  (4.6  and  4.8)  for  the  induced  circulation.  These 
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terms  signify  that  the  mean  acceleration  of  the  induced  circulation, 
in  general,  is  non-zero  and  is  influenced  by  the  detailed  time  history 
of  the  fluctuations  of  the  eddy  flux  convergences  (divergences)  and 
not  merely  by  their  mean  values.  The  circumstances  under  which  the 
acceleration  terms  will  go  to  zero  is,  in  fact.,  the  content  of  the 
so  called  non-acceleration  theorem  discussed  by  several  authors  in 
slightly  different  contexts  (Andrews,  1985;  Charney  and  Drazin,  1961; 
Dunkerton,  1978,  1980.) 
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FIGURE  CAPTIONS  FOR  CHAPTER  IV 


Fig.  4.1.  The  mean  linear  circulation  (see  text  for  definition). 

Fig.  4.2.  The  mean  induced  circulation  (see  text  for  definition). 

Fig.  4.3.  Heat  sources  and  sinks  due  to  the  mean  Eulerian  circula¬ 
tion  expressed  as  equivalent  temperature  changes.  Notice  that  the 
source  regions  have  been  eliminated  from  this  figure  and  from  Figures 
4.4,  4.5,  and  4.6. 

Fig.  4.4.  Same  as  Figure  4.3  except  due  to  the  mean  linear  circu¬ 
lation  . 

Fig.  4.5.  Same  as  Figure  4.3  except  due  to  the  mean  induced  circu¬ 
lation  . 

Fig.  4.6.  Same  as  Figure  4.3  except  due  to  the  mean  eddy  heat  flux 
convergences  (divergences) . 
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LINERR  WINDS  12  HR  RVE 
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TEMP  CHANGES  EDDY  HEAT  FLUX 
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CHAPTER  V 


LAGRANGIAN  MOTIONS,  EULERIAN  MOTIONS,  AND  STOKES  DRIFT 

5.1  TRAJECTORY  CALCULATIONS 

As  described  in  Chapter  3,  the  storm-time  winds  vary  both  in  space 
and  time.  Because  of  these  variations,  the  actual  trajectory  of  an 
identified  fluid  marker,  referred  to  as  the  Lagrangian  trajectory, 
will  be  different  from  the  mean  Eulerian  trajectory  along  the  stream¬ 
lines  corresponding  to  the  time-averaged  wind  vectors  at  f  xed  points 
in  space. 

In  this  paper,  the  Lagrangian  trajectory  calculations  are  perform¬ 
ed  using  an  iterative  method.  After  fixing  the  starting  position  of 
the  fluid  marker  on  our  model  grid,  we  use  the  initial  velocity  field, 
V ( t Q ) ,  at  that  position  to  calculate  the  displacement,  V(tQ)At,  where 
At  is  a  time  increment  of  10  minutes.  We  next  find  the  new  grid 
location  of  the  marker  at  tQ  +  At,  where  its  velocity  is  V(t0  +  At). 
V ( t q )  and  V ( t q  +  At)  are  then  averaged  and  used  to  obtain  a  better 
value  for  the  displacement  of  the  marker  from  its  initial  position. 
This  procedure  is  repeated  until  the  displacement  from  iteration  to 
iteration  does  not  vary  significantly.  We  then  use  the  marker* s  new 
position  as  the  starting  point  for  calculating  the  displacement  for 
the  next  time  increment.  The  trajectories  are  calculated  in  this 
manner  for  a  maximum  time  period  of  12  hours,  and  several  such  trajec¬ 
tories  are  shown  in  Figures  5.1  and  5.3. 

Figure  5.1  shows  the  Lagrangian  trajectory  of  the  marker  started 
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at  400  km  and  70*  colatitude.  In  this  example  it  is  assumed  that  the 
thermosphere  is  at  rest  except  for  the  motion  associated  with  the 
storm-generated  waves.  Five  such  waves  are  shown  to  influence  the 
motion  of  the  marker  in  the  12-hr  period  of  storm  simulation.  At  the 
end  of  this  period,  the  marker  is  found  to  have  moved  a  horizontal 
distance  of  approximately  1.5  degrees  in  colatitude  (-150  km)  and 
upward  by  about  3  km.  This  displacement  of  the  marker  between  its 
starting  and  ending  positions  is  marked  by  an  arrow,  labelled  L,  and 
is  its  Lagrangian  displacement. 

Figure  5.2  shows  the  Lagrangian  trajectories  of  the  marker  started 
at  a  different  location  (110*  colatitude,  257  km  altitude).  In  this 
case,  the  marker  is  influenced  both  by  the  waves  and  the  storm-time 
meridional  circulation.  We  can,  therefore,  distinguish  between  the 
net  Lagrangian  displacement  and  the  net  Eulerian  displacement  of 
the  marker  in  a  given  time  interval.  In  Figure  5.2a  this  time 
interval  is  chosen  to  be  the  time  of  passage  of  the  first  wave,  and 
in  Figure  5.2b,  of  the  first  two  waves.  In  each  case,  the  net 
Lagrangian  and  net  Eulerian  displacements  are  indicated  by  the 
vectors  L  and  E,  repectively. 

Figure  5.3  gives  examples  of  the  Lagrangian  and  Eulerian  motions 
of  markers  started  at  the  same  altitude  (-257  km)  but  at  different  co¬ 
latitudes  (40  * ,  80 * ,  110 * ) .  Figure  5.3a  is  for  the  case  for  storm-gene¬ 
rated  waves  and  circulation  alone,  whereas  Figure  5.3b  includes  the 
seasonal  background  meridional  circulation  as  well.  Figure  5.3a  illu¬ 
strates  the  following  features:  1)  Individual  air  parcels  originat- 
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ing  near  the  heat  source  region  do  not  have  time  to  reach  the  equator 
in  the  12-hr  duration  of  the  storm*  They  can  at  most  travel  a  horizon¬ 
tal  distance  of  15'  in  colatitude  (-1500  km)  .  2)  At  intermediate 
colatitudes,  air  parcels  move  generally  downward  and  equatorward. 

3)  Near  the  equator,  the  air  parcel  transport  is  predominately  down¬ 
ward,  and  is  initiated  with  a  time  delay  of  approximately  2  to  3  hours 
after  the  beginning  of  the  storm,  signalling  the  first  arrival  of  the 
storm-generated  waves. 

The  differences  in  the  above  description,  brought  about  by  the 
inclusion  of  the  seasonal  background  winds,  are  the  larger  meridional 
displacements  (-35*  in  colatitude  near  the  heat  source  region)  of 
air  parcels,  and  the  shift  in  location,  away  from  the  equator,  of  the 
region  of  maximum  vertical  descent. 

5.2  STOKES  RELATIONS 

Following  Stokes  (1847)  and  Longuet-Higgins  (1969),  we  can  relate 
the  difference  between  the  net  Lagrangian  and  the  net  Eulerian 
displacements  to  the  Stokes  drift  (Figure  5.2).  To  obtain  this 
relation,  we  start  with: 

-4-4  -4  *4 

V  (  r,  t  )  -  V  (  r0,  t)  +  V  (  r0,  t  )  *  V  v  (  r0,  t  )  dt  (5.1) 

Jt0 

-4  -4  ->  -4 

where  V(  r,t)  is  the  Lagrangian  velocity,  V(  rQ,t)  is  the  Eulerian 

velocity  and  the  difference  between  them  is  the  Stokes  drift  velocity 

-4  -4 

appearing  under  the  intergral  sign.  Here  rg,  r  are  the  position 


88 


vectors  of  the  fluid  marker  during  a  time  interval  between  tQ  and  t 

~4 

with  A  r  =  r  -  tq  ,  If  we  assume  that  A  r  is  small  compared  to  the 
local  spatial  scale  of  the  velocity  variations,  we  can  write: 


-4  -4  -4  -4  -4  -4  *4 

V  (  r,  t)  =  V  (  rg,  t)  +  A  r  •  V  V  (  rg,  t) 


A  r  =  V  (  rg,  t  )  dt 
Jt0 


from  which  equation  5.1  follows.  To  obtain  the  Stokes  velocity 
relation,  we  take  the  mean  over  the  period  (t-tg),  denoted  by  an 
overbar,  and  write: 


-4  4 

V  <  r,  t  ) 


*4  4  ft  4  4 

-  V  (  rg,  t  )  +  V  (  rg,  t  ) 


4  4 

V  V  (  r0/t  >  dt  (5.4) 


-4  -4  4 

VT  =  V  + 

L  E  S 


-4  “4  -4 

where  V  ,  V  ,  V  are  the  mean  Lagrangian  velocity,  the  mean  Eulerian 
L  E  S 

velocity,  and  Stokes  drift  velocity,  respectively.  In  Figure  5.2,  we 

■4  «  -4  *4  p  “4 

have  shown  the  net  displacements  with  L  *  J  dt,  E  *  J  V£  dt,  and 

•4  p  4 

S  =  dt. 

o 

Our  trajectory  calculations  can  be  used  to  account  for  the  earli¬ 
est  temperature  changes  observed  near  the  equator  about  two  to  three 
hours  after  the  beginning  of  the  storm  (figure  3.2b).  These  tempera¬ 
ture  changes  could  not  have  been  due  to  the  material  transport  of 


fluid  fron  the  heat  source  regions.  On  the  other  hand,  the  earliest 


equatorial  temperature  signal  accompanies  the  arrival  of  the  first 
storm-generated  wave.  This  wave  runs  ahead  of  the  meridional  circu¬ 
lation  that  gets  established  behind.  Both  the  waves  and  the  circu¬ 
lation  are  agents  responsible  for  the  descent  of  air  parcels  and  the 
consequent  heating  of  the  equatorial  thermosphere  during  the  storm. 
This  conclusion  is  consistent  with  the  role  of  the  induced  circulation 
in  the  equatorial  region,  discussed  in  Chapter  4. 
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FIGURE  CAPTIONS  FOR  CHAPTER  V 


Fig.  5.1.  Trajectory  of  fluid  marker  solely  due  to  storm-generated 
wave  disturbances.  The  marker  is  started  at  70*  colatitude  and  an 
altitude  of  400  km.  The  numbers  1  to  5  correspond  to  the  beginning 
times  of  the  five  disturbances  influencing  the  fluid  marker.  Triangles 
show  the  locations  of  the  marker  at  2-hour  intervals.  The  vector  L 
denotes  the  net  Lagrangian  cr  ^placement . 

Fig.  5.2.  Trajectory  of  the  fluid  marker  under  the  influence  of  storm¬ 
time  waves  and  meridional  circulation  started  at  257  km  and  110*  cola¬ 
titude.  a)  The  Lagrangian  trajectory  (solid  line)  is  shown  during  the 
passage  of  the  first  storm-generated  wave  with  the  locations  of  the 
marker  indicated  by  triangles  at  30  minute  intervals.  In  addition, 
the  hypothetical  positions  of  the  marker,  if  it  had  moved  with  the 
mean  Eulerian  velocity  during  this  interval,  are  also  indicated  with 
asterisks  at  30  minute  intervals.  The  vectors  L,  E,  and  S  are  for 
the  net  Lagrangian  displacement,  the  net  Eulerian  displacement,  and 
the  Stokes  drift,  respectively.  b)  The  same  trajectories  and  vectors 
are  shown  during  the  passage  of  the  first  two  waves.  Note  that  in 
this  case  the  mean  Eulerian  velocity  is  obtained  by  averaging  over 
two  wave  periods . 

Fig.  5.3.  Examples  of  Lagrangian  and  Eulerian  motions  of  markers 
started  at  257  km  and  at  colatitudes  of  40*,  80*,  and  110*.  The  solid 
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lines  are  the  Lagrangian  trajectories  and  the  dashed  lines  are  the 
Eulerian  trajectories.  The  asterisks  and  triangles  in  these  trajec¬ 
tories  are  at  2-hr  intervals.  The  triangles  are  labelled  for  storm- 
time  in  hours  to  aid  in  following  the  marker's  path.  a)  For  the  case 
with  storm-generated  waves  and  circulation  alone.  b)  For  the  case 
that  also  includes  the  seasonal  background  meridional  circulation. 
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CHAPTER  VI 


DISCUSSION  AND  CONCLUSIONS 


In  this  paper,  we  have  studied  the  waves  and  the  global  meridional 
circulation  generated  by  a  geomagnetic  storm.  In  particular,  we  have 
investigated  the  interaction  between  the  waves  and  the  mean  flow  using 
a  two-dimensional  pole-to-pole  model  in  which  the  solstitial  background 
winds  have  been  included. 

The  model  storm  that  we  have  used  consists  of  a  series  of  substorms 
which  produced  four  distinct  heat  pulses  lasting  through  the  12-hour 
period  of  storm  simulation.  Although  we  have  specific  examples  of  such 
storms,  as  shown  in  Figure  2.5a,  this  type  of  storm  associated  with 
sustained  wave  activity  may  not  be  typical.  In  another  well-documented 
class  of  storms,  large  amplitude  waves  are  generated  mainly  in  the 
sudden  commencement  phase  {Roble  et.  al.,  1978). 

In  our  model  storm,  we  are  able  to  identify  six  distinct  wave  dis¬ 
turbances  corresponding  to  the  four  pulses  in  energy  input.  These 
wave  disturbances  are  observed  to  consist  of  a  spectrum  of  waves  which 
travel  equatorward  from  their  sources  of  origin.  When  the  disturbances 
are  observed  to  be  propagating  poleward,  their  origin  can  be  traced  to 
the  opposite  hemisphere. 

The  amplitudes  of  the  disturbances  in  middle  latitudes  range  be¬ 
tween  50-100  m/s  in  meridional  velocity  and  have  a  maximum  of  7.5  m/s 
in  vertical  velocity.  Their  phase  speeds  can  be  as  high  as  900  m/s. 
The  temperature  fluctuations  associated  with  the  disturbances  are  on 
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the  order  of  50*  K.  These  wave  characteristics  are  similar  to  what 
have  been  reported  from  observations  (Hernandez  and  Roble,  1978;  Rees 
et .  al.,  1984;  Roble  et .  al.,  1978). 

The  apparent  periods  of  the  disturbances  are  found  to  increase 
with  distance  from  the  wave  source,  indicating  that  the  disturbances 
are  dispersive.  The  disturbance  amplitudes  are  reduced  and  the  phase 
speeds  are  Doppler-shifted  in  the  direction  of  the  background  wind. 
However,  the  phase  shifts  as  well  as  the  amplitude  reductions  are 
small  due  to  the  fact  that  the  phase  speeds  of  the  disturbances  (450- 
900  m/s)  are  large  compared  to  the  background  wind  speeds  (10-100  m/s)  . 
The  phase  speeds  of  the  disturbances  are  also  influenced  by  the  tempe¬ 
rature 

The  disturbances  tilt  with  height  in  a  manner  consistent  with  the 
propagation  of  energy  away  from  the  high-latitude  heat  source  and  up¬ 
ward.  They  are  asymmetric  in  the  sense  that  over  a  period  of  the  dis¬ 
turbance,  the  average  velocity  or  net  temperature  change  is  not  zero. 
Moreover,  the  phase  relationships  between  the  disturbance  parameters 
are  such  that  over  the  period  of  the  storm,  the  net  transport  of  both 
meridional  momentum  and  heat  is  equatorward. 

The  results  of  the  interaction  between  the  waves  and  the  mean  flow 
can  be  portrayed  in  terms  on  the  induced  circulation.  Our  results  for 
the  induced  circulation  show  that  its  influence  is  clearly  detectable 
in  the  time-averaged  equivalent  temperature  changes  due  to  the  storm. 
In  the  equatorial  region  at  high  altitudes  these  changes  are  about  1/4 
to  1/3  of  the  changes  due  to  the  mean  Eulerian  circulation. 
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The  spatial  and  temporal  variations  of  the  storm-time  winds  are 
responsible  for  the  difference  between  the  Lagrangian  and  mean  Eulerian 
trajectories  rn  the  storm-time  thermosphere.  Lagrangian  trajectory 
calculations  clearly  show  that  there  cannot  be  material  transport  of 
fluid  from  the  source  latitudes  to  the  equator.  The  earliest 
equatorial  motions  are  induced  by  storm-generated  waves  reaching  the 
equator. 
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